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ABSTRACT 

The WFC3 Infrared Spectroscopic Parallel Survey (WISP) uses the Hubble Space Telescope infrared 
grism capabilities to obtain slitless spectra of thousands of galaxies over a wide redshift range including 
the peak of star formation history of the Universe. We select a population of very strong emission- 
line galaxies with rest-frame equivalent widths higher than 200 A. A total of 176 objects are found 
over the redshift range 0.35 < z < 2.3 in the 180 arcmin 2 area that we have analyzed so far. This 
population consists of young and low-mass starbursts with high specific star formation rates (sSFR). 
After spectroscopic follow-up of one of these galaxies with Keck/LRIS, we report the detection at 
z = 0.7 of an extremely metal-poor galaxy with 12+Log(0/H)= 7.47 ± 0.11. After estimating the 
AGN fraction in the sample, we show that the high-EW galaxies have higher sSFR than normal 
star-forming galaxies at any redshift. We find that the nebular emission-lines can substantially affect 
the total broadband flux density with a median brightening of 0.3 mag, with some examples of line 
contamination producing brightening of up to 1 mag. We show that the presence of strong emission 
lines in low-z galaxies can mimic the color-selection criteria used in the z ~ 8 dropout surveys. In order 
to effectively remove low redshift interlopers, deep optical imaging is needed, at least one magnitude 
deeper than the bands in which the objects are detected. Without deep optical data, most of the 
interlopers cannot be ruled out in the wide shallow Hubble Space Telescope imaging surveys. Finally, 
we empirically demonstrate that strong nebular lines can lead to an overestimation of the mass and 
the age of galaxies derived from fitting of their SED. Without removing emission lines, the age and 
the stellar mass estimates are overestimated by a factor of 2 on average and up to a factor of 10 for 
the high-EW galaxies. Therefore the contribution of emission lines should be systematically taken 
into account in SED fitting of star-forming galaxies at all redshifts. 

Subject headings: galaxies: evolution - galaxies: statistics - galaxies: high-redshift - infrared: galaxies 
- surveys - cosmology: observations 



1. INTRODUCTION 

The characterization of galaxies undergoing their first 
major star formation episode is of great interest in mod- 
ern astrophysics. This class of galaxies is expected to 
host young and massive stars ionizing the interstellar 
medium, and therefore to exhibit strong nebular emis- 
sion lines. However, many efforts dedicated to the un- 
derstanding of galaxy formation and evolution are in gen- 
eral limited by different factors or colored by the way in 
which galaxy samples are assembled. Broadband sur- 
veys yield magnitude-limited samples that are biased 
towards bright contin uum objects such as the Lyman 
Break Galaxies (LBGs. St eidel et al.lll996llShaplev et al.1 
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20031 iVanzella et al.l2009tlHathi et all2010l) , and can fa- 
vor relatively massive galaxy populations. On the other 
hand, ground-based searches for rest-frame optical emis- 
sion lines from high-redshift galaxies are severely im- 
pacted by the bright NIR background. Observing in 
small wavelength windows between OH terrestrial air- 
glow will significantly restrict the survey volume. 

In this context, space-based observations offer consid- 
erable advantages for the detection of z ~ 1 — 2 emission- 
line galaxies. The Wide field Camera 3 (WFC3) on- 
board the Hubble Space Telescope (HST) has dramati- 
cally improved the power of slitless spectroscopy from 
above the atmosphere. Beginning in cycle 17, we are con- 
ducting the WFC3 Infrared Spectroscopic Parallel Sur- 
vey (WISP). In about 500 orbits of parallel observations, 
the program will obtain slitless spectra of galaxies at the 
peak of the st ar formation history of the Universe (e.g. 
IHo D kinsl[200l . The G 102 and G 14 i LR grisms offer con- 
tinuous wavelength coverage from 0.8 to 1.7 /xm, allowing 
the selection of emission-line galaxies at 0.35 < z < 2.5, 
and potentially bright Lycc emitters at z > 6.5 (see 
lAtek et all 120101 for a complete description of the sur- 
vey). 

In this paper, we present a sample of emission-line ob- 
jects with extremely high equivalent widths. The equiv- 
alent width of recombination lines indicates the ratio of 
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the current star formation rate (SFR) to t he past aver- 
age SFR of a galaxy (e.g. iKennicu ttl TT998h . Therefore, 
the high-EW selection picks up galaxies with a strong 
ongoing star formation episode. 

Previously, ground-based narrow-band surveys were 
used to identify a large s ample of ultra strong emis- 
sion l ines galaxies (USELs. lKakazu et al.ll2007t IHu et al.l 
2009), but were restricted to a narrow range of redshifts 
up to z = 1. In the local universe, iCardamone et al.l 
(2009) found in the Sloan Digital Sky Survey (SDSS) 
compact star-forming galaxies, called "green peas" , that 
have very blue r — i colors because of their very strong 
[Oin]A5007 emission line. Thanks to the unprecedented 
sensitivity of the WFC3/IR detectors, the depth of the 
direct imaging combined with the resolution and the 
wavelength coverage of the grism spectroscopy allows us 
to uncover this particular class of objects independently 
of the continuum brightness over a wide redshift range. 

The selection of high-EW star-forming galaxies with 
the WFC3 grisms probes a lower mass range than previ- 
ous studies at z > 1 (e.g. lErb et al.1 12006a iDaddi et al.l 
120071 ) . This allows us to compare the star-formation 
efficiency of these dwarf galaxies with normal star- 
forming galaxies at z ~ 1 — 2. Indeed, several stud- 
ies have established a correlation between the star 
formation rate and the stellar mass up to redshift 
z = 7 (e.g. iBrinchmann fc Ellidl2000t lElbaz et alj|2007t 
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20091 lOliver et alJl2010t lLabbe et al.l 120101 ) . The shape 



and evolution of this relationship is interpreted as 
the result of smo oth gas accretion (|Bouche et al.l 1201a 
iDave et aT1l201lD . However, at high redshifts, the stud- 
ies focused on a specific mass range due to the difficulty 
of selecting faint-continuum galaxies. 

In a broader context, we also discuss the implications 
of these strong line-emitters on the selection of z ~ 8 
galaxies and on the physical properties of galaxies de- 
rived from SED modeling. 

We present in Section [2] the WISP survey observa- 
tions and the data reduction. Section [3] is devoted to 
the follow-up observations. The selection of the high- 
EW population is presented in Section |H In Section |6] 
we discuss the mass and star-formation properties of the 
high-EW galaxies. We present in Section [7] the metallic- 
ity measurement of one example and put it into context. 
We investigate the contribution of nebular lines to the 
total broadband flux density and their astrophysical im- 
plications in Sectionl5land [5~2l Our conclusions are given 
in Section [5] Throughout, we assume a A-dominated flat 
universe, with Hq = 70 km s _1 Mpc -1 , fl\ — 0.7, and 
ft m = 0.3. All magnitudes are in AB system. 

2. THE WISP SURVEY 

The WISP survey exploits the powerful pure paral- 
lel mode of the HST to obtain slitless spectra of nearly 
200 uncorrelated fields in 500 orbits within two observing 
programs (PI = M. Malkan): GO 11696 and 12283 (In- 
frared Survey of Star Formation Across Cosmic Time). 
During long int egrations with the Co smic Origin Spec- 
trograph (COS. fOsterman et al.ll201ll ) o r the Space Tele- 
scope Imager and Spectrograph (STIS, [Woodgat e et al 



1998), the WISP program uses WFC3 (|Kimble et al 



20081) to observe an adjacent field in parallel, at a fixed 
offset of about 5' from the primary target. The data in- 



clude slitless spectroscopy with G102 (0.8 - 1.2 fiia) and 
G141 (1.1-1.7 fim) grisms, and direct NIR imaging with 
the F110W and F140W filters. As of August 2010, we 
added the F475X and F600LP filters of the WFC3/UVIS 
channel to our observing program, with typical exposure 
times of 400 sec in each band. The fields analyzed in 
this paper do not have UVIS observations. For the long 
visits (typically 4 orbits or more), we observe with both 
grisms with a typical exposure time ratio of G102:G141 
= 2.5:1. For shorter visits the fields are observed in G141 
only. The direct images are obtained in the same orbits 
as the grism imaging with a ~ 6:1 grism:direct integra- 
tion time ratio. It is necessary for the spectral extraction 
from slitless data to use the direct images to provide the 
input catalog of objects. The object positions, sizes and 
shapes are then used, in conjunction with a dispersion 
solution that translates the position from the direct im- 
age to the spectral trace in the grism one, to extract flux 
and wavelength calibrated 2D and ID spectra. 

In this paper, we analyzed a total of 54 fields presented 
in Table [T] among which two have optical follow-up (cf. 
Section [3]). We note that our effective area represents 
about 70% of the full WFC3/IR frame because of two 
things. First, objects that fall outside the right edge of 
the direct image frame will have their zero order spectra 
falling in the right side of the grism frame (zeroth order 
is dispersed to the left side of the object position). Thus, 
we cannot correct the zeroth order contamination in that 
part of the detector. Second, the left edge of the grism 
frame will receive spectra from objects that are outside 
the direct image (first order dispersed to the right side of 
the object position), preventing any flux or wavelength 
calibration for those sources. 

All the data were processed with the WFC3 pipeline 
CALWF3 (version 2.1) to correct for bias, dark, flatficld 
and gain vari ations. Then, the sl itless extraction pack- 
age aXe 2.0 (|Kiimmel et al.l I2009I) is used for the spec- 
tral extraction. Spectra with important contamination 
are removed from the analysis because the contamina- 
tion estimate of aXe is not sufficiently accurate. A com- 
pl ete descript i on of the data reduction steps is presented 
in I Atek et all (|2010f ). 

3. OPTICAL FOLLOW-UP SPECTROSCOPY AND IMAGING 

We have obtained optical spectra of objects in two 
adjacent WISP field s with the Low Resolution Imaging 
Spectrometer rLRIS. lQfeltllllTQMISteidel et al.ll200"l 
on the Keck I telescope. The two WFC3 fields, WISP5 
and WISP7 in Table [TJ are immediately adjacent to each 
other and fit within one LRIS slit mask. We assigned 
1.2" wide slits to 13 emission line objects, 5 of which be- 
long to the high-EW sample. Three exposures of 1800s 
each were used for a total integration of 5400s. The 400 
1/mm grism blazed at 3400 A was used on the blue arm 
and the 600 1/mm grating blazed at 1.0 fim was used on 
the red arm for a pixel scale of ~ 1.1 A pix -1 and ~ 0.8 A 
pix -1 , respectively. The plate scale is 0.135" pix -1 and 
the spectral resolution for an object that fills the slit is 
about 8.1 A and 5.6 A in the blue and red side, respec- 
tively. The seeing during observations was around 1". 
The spectrophotometric standard Wolf 1346 was used to 
flux-calibrate the spectra. 

We also observed these two fields on April 20 2010 
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Fig. 1. — Strong emission-line galaxies across a wide range of redshift in the WISP Survey. Left: Montage spectra of 176 high-EW 
emission-line galaxies stacked in the Y-axis direction with increasing redshift. Consequently, we see the emission lines we routinely detect 
in WISPS being shifted in the wavelength dimension. We note that the slope of the emission-line curve is less steep in some parts of this 
plot. i.e. we detect less objects per redshift bin, which is due the decreasing sensitivity towards the edge of the grisms. Right: An example 
of a galaxy at z ~ 0.7 showing strong emission lines. The direct F110W image and the G102 and G141 2D spectra are shown above the 
calibrated ID spectrum. The J and H magnitudes are also given in the inset. The object shows bright Ha and [Oin] lines at 0.85 and 1.12 
lira, respectively, on top of a very faint continuum. The He* line is detected in both grisms, in the overlapping region. 



to obtain optical images with the Large Format Camera 
(LFC) on the Hale 5m telescope at Palomar observatory. 
The two fields were easily covered by the mosaic camera 
with a field of view of 24 arcmin diameter. We used the g' 
and i' filters with a total exposure time of 1500 and 3000 
sec, respectively. On readout, we binned 2x2, resulting 
in a pixel scale of 0.36" per pixel. The typical seeing 
during observations was about 1". We adopted a 5 point 
dithering pattern with a 15" displacement. A summary 
of the observations is given in Table [T] 

The images were corrected for bias and flatfield using 
standard IRAF reduction packages. The flux calibration 
was performed by comparing the photometry of the field 
stars of our im ages with the thos e of the Sloan Digital Sky 
Survey fSDSS lYork et alJl2000T ) that we obtained from 
the data release DR8. Finally, we aligned the ground- 
based frames to our IR ones using GE0MAP and GE0TRAN 
tasks in IRAF. In order to measure the object magnitudes 
in the optical images we need to match the photome- 
try aperture between ground-based and the HST images. 
We first convolve the WFC3 images with a gaussian ker- 
nel to match the PSF of the Palomar data. Then, we 
measure the flux in the convolved images using an aper- 
ture of 5 WFC3/IR pixels. The aperture size adopted is 
large enough to encompass most of the object flux for our 
compact sources, without being contaminated by nearby 
sources. We then compare the measured magnitudes to 
the AUTO magnitudes measured with SExtractor in the 
original WFC3 images. This gives us the aperture correc- 
tion to apply to the magnitudes measured in the optical 
data. 

4. HIGH EQUIVALENT- WIDTH POPULATION 

The WISP survey is particularly sensitive to high spe- 
cific SFR and strong emission-line galaxies. In the fields 
we observed so far, we often see isolated emission lines 
in the grism frames, with either a faint continuum or no 
continuum at all. Such objects are easily identifiable by 
visually inspecting the final ID spectra. We only search 



part of the field, for which the zeroth order information 
is available. To this aim, we have developed a flagging 
scheme in our reduction pipeline that marks the position 
of the zeroth order in the 2D and ID spectra, to prevent 
their identification as emission lines. We then use a cus- 
tom IDL procedure based on a least square method to 
fit all the lines in each spectrum with a Gaussian and a 
polynomial continuum model. 

After the line measurement we have retained only 
galaxies with rest-frame EW higher than 200 A. Exam- 
ples of strong emission-line objects in our catalog are 
presented in Figure [T] The strong lines are typically 
Ha and [Om]. The sample reaches a 3<r flux limit of 
7 x 10 -17 erg s _1 cm" 2 over a 180 arc min 2 area, which i s 
about three times the area surveyed in lAtek et al.l {2010). 
Both IR grisms are used to observe 24 fields (80 arcmin 2 ) 
over a wavelength range of 0.8-1.7 /jm, and 30 fields (100 
arcmin 2 ) are observed with G141 only covering a wave- 
length range 1.1-1.7 /im. In Figure [2] we present the 
rest-frame equivalent width distribution of the high-EW 
sample. It is well known that the Ha equ ivalent width is 
an in dicator of the age of a galaxy (e.g. iLeitherer et al.l 
1999). The Ha equivalent width indicates the ratio of the 
current SFR to the lifetime averaged SFR, also known a s 
the birthrate parameter b (|Scalolll986t lKennicutt]|1998lL 
The continuum luminosity is representative of the low 
mass stars with a lifetime of several Gyrs. The strength 
of Ha line evolves on much shorter timescales, as it is the 
result of hot, short lived O and B stars. A galaxy with a 
constant star formation will have a birthrate parameter 
of unity, while b > 1 indicates that the current episode of 
star formation is stro nger than the averag e SFR. Using 
Starburst99 models (jLeitherer et al.l["l999t ). we find that 
the 200 A limit on EW(Ha) selects either young galax- 
ies (less than 10 Myr) with instantaneous burst episode, 
or galaxies with continuous star formation older than a 
burst and less than 1 Gyr. 

We note that for extremely high-EW lines, the con- 
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TABLE l 

Summary of WFC3 Imaging and Spectroscopy 



Field 


RA 


DEC 


F110W 


G102 


F140W 


G141 




(HMS) 


(DMS) 


(sec) 


(sec) 


(sec) 


(sec) 



WISP1 


01 


06 


35.29 


+15 


08 


53.8 


884 


4815 


506 


2609 


WISP2 


01 


25 


10.02 


+21 


39 


13.7 








328 


1906 


WISP5* 


14 


27 


06.64 


+57 


51 


36.2 


1034 


5515 


1034 


5515 


WISP6 


01 


50 


17.18 


+13 


04 


12.8 


609 


3609 


862 


5015 


WISP7* 


14 


27 


05.91 


+57 


53 


33.7 


834 


6318 


1112 


6224 


WISP8 


11 


51 


51.62 


+54 


40 


41.2 


1662 


9021 


581 


2509 


WISP9 


12 


29 


44.31 


+07 


48 


23.5 


759 


4612 


684 


3712 


WISP10 


09 


25 


07.84 


+48 


57 


03.0 


631 


3909 


406 


2209 


WISP11 


11 


02 


17.38 


+10 


54 


25.4 


556 


3709 


456 


2006 


WISP12 


12 


09 


25.25 


+45 


43 


19.8 


1312 


8221 


606 


3009 


WISP13 


01 


06 


38.77 


+15 


08 


26.2 


556 


3009 


506 


2409 


WISP14 


02 


34 


56.80 


-04 


06 


54.5 


834 


6215 


181 


2809 


WISP15 


11 


09 


42.47 


+26 


21 


56.0 


1612 


8321 


531 


2609 


WISP16 


02 


34 


54.72 


-04 


06 


42.5 


1087 


6921 


584 


2509 


WISP18 


12 


29 


17.25 


+10 


14 


00.6 


534 


3512 


762 


3824 


WISP19 


02 


34 


54.29 


-04 


06 


30.5 


1187 


8721 


484 


2809 


WISP20 


14 


09 


41.15 


+26 


22 


15.1 


1815 


8430 


559 


2812 


WISP21 


09 


27 


55.77 


+60 


27 


05.3 








353 


2006 


WISP22 


08 


52 


44.99 


+03 


09 


09.6 








253 


1806 


WISP23 


09 


43 


16.12 


+05 


27 


37.1 








681 


4115 


WISP24 


12 


18 


41.90 


+29 


■52 


51.8 








278 


1806 


WISP25 


10 


08 


42.49 


+07 


11 


10.3 








734 


4115 


WISP26 


08 


45 


17.91 


+22 


54 


58.5 


840 


5521 


384 


2209 


WISP27 


11 


33 


05.98 


+03 


28 


02.2 


1040 


6218 


384 


2206 


WISP28 


09 


35 


46.24 


+14 


27 


48.5 








634 


3515 


WISP29 


12 


02 


59.55 


+48 


05 


21.2 








328 


1906 


WISP30 


10 


28 


18.91 


+39 


17 


14.1 








253 


1806 


WISP31 


08 


43 


27.99 


+26 


16 


39.7 








659 


4412 


WISP32 


13 


05 


20.67 


-25 


38 


05.5 








181 


3512 


WISP33 


10 


00 


02.78 


+12 


14 


56.4 


959 


6415 


484 


2809 


WISP34 


09 


34 


58.06 


+02 


02 


03.7 








709 


4215 


WISP35 


13 


03 


46.90 


+29 


53 


03.8 








634 


3812 


WISP36 


13 


10 


31.52 


+41 


23 


10.3 


1237 


8621 


584 


2809 


WISP38 


12 


25 


13.33 


-02 


49 


08.4 








709 


4312 


WISP39 


10 


09 


38.23 


+30 


00 


48.7 








328 


1706 


WISP40 


02 


16 


20.12 


-39 


02 


28.1 








584 


4112 


WISP41 


12 


08 


28.31 


+45 


38 


48.2 


1084 


6615 


584 


2909 


WISP42 


10 


01 


04.36 


+50 


21 


02.8 


1812 


10424 


784 


4012 


WISP43 


21 


04 


07.62 


-07 


23 


00.8 


909 


6315 


434 


2809 


WISP44 


11 


12 


15.18 


+35 


36 


55.5 








734 


4515 


WISP45 


12 


37 


26.68 


+01 


25 


18.2 








796 


4812 


WISP46 


22 


37 


58.48 


-18 


42 


05.4 








328 


1906 


WISP46 


22 


37 


58.48 


-18 


12 


05.4 








328 


1906 


WISP47 


13 


19 


32.40 


+27 


26 


40.0 








587 


3212 


WISP49 


14 


44 


44.84 


+34 


27 


45.9 


1093 


6018 


406 


2406 


WISP50 


22 


22 


20.25 


+09 


36 


36.8 








328 


1906 


WISP51 


15 


13 


13.77 


+36 


33 


22.0 








278 


2006 


WISP52 


13 


30 


22.62 


+28 


11 


01.0 








634 


4218 


WISP53 


15 


11 


12.91 


+40 


25 


39.3 








253 


1806 


WISP54 


15 


44 


57.26 


+48 


15 


23.2 








303 


2006 


WISP56 


16 


16 


50.58 


+06 


36 


43.6 








684 


4312 


WISP57 


12 


33 


16.74 


+47 


53 


27.5 








659 


4615 


WISP59 


15 


50 


22.59 


+39 


59 


19.2 








734 


4515 


WISP61 


23 


09 


00.35 


-09 


07 


17.7 








709 


4312 



Note. — Observation information for the fields analyzed in this paper. WISP5 and 
WISP7 have optical imaging in g' and i' bands with exposure times of 1500 and 3000 sec, 
respectively. The two fields were also covered by a slitmask with LRIS at Keck I for a total 
exposure time of 5400 sec. 1.2"wide slits were used. 



tinuum remains undetected in the grism spectra, mak- 
ing the EW measurement highly uncertain. Given the 
low spectral resolution of the WFC3 grism we are not 
able to separate the [Nil] lines from Ha. However, be- 
cause the high-EW galaxies are likely metal-deficient (see 
Section [7]) , we do not expect a large contribution from 
the [Nil] lines. In fact, given an [Olll] A5007/H/3 ra- 
tio, one can infer the [ Nil] /Ha ratio from t he [Sii]/Ha 
one. From Figure 1 of lKewlev et~all (f2006h . and using 
the measured [Sn]/([Nn]+Ha), we derive the [Sii]/Ha 



and then [Nil] /Ha for 17 high-EW galaxies. We found a 
median contribution of [Nil] to Ha of about 8 %. 

The strong [Oni] A5007 line can also be the result of 
a photoionization dominated by harder radiation field 
from an Active Galactic Nucleus (AGN). This will fa- 
vor collisionally-excited lines rather than recombination 
lines. In order to estimate the AGN contribution in the 
sample we use the BPT diag nostic diagram proposed by 
iVeilleux fc Osterbrockl (fl987l) using [Olll] A5007/H/3 ver- 
sus [Sli] A(6717+6732)/Ha line ratios, where Ha is cor- 
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rected for [Nil] contamination following the procedure 
described above. We note that using the BPT diagram 
to classify low-metallicity dwar f galaxies with high sSFR 
can be somewhat ambi guous (iBrinchmann et al.l 120081 : 
llzotov fe Thuanll2008D . IBrinchmann et al.l (|2008T) show 
that high-sSFR galaxies may have higher values of the 
ionization parameter than normal star-forming galaxies, 
and tend to be offset in the BPT diagram towards the 
AGN region. 

The classification of the six objects is presented in Fig- 
ured] We identified one AGN in a subsample of six galax- 
ies that fall in the right redshift range (0.7 < z < 1.4) 
and have detectable H/3 and [Sn] lines. This represents 
a fraction of about 17 %, although the statistics are 
too small to derive a secure AGN fraction for the en- 
tire sample. This is cons istent with the fraction found in 
z ~ - 3 — 2 Lyq emitters (lAtek et alJl2009HScarlata et al.1 
[200l iNilsson et al J 1200a iCowie et all 12011 1201 lh . and 
what is locally fo und in emission-line galaxy samples (e.g. 
IHao et al.l[2005h . 

This high-EW selection is co mparable to that o f 
ground-based narrowband surveys dKakazu et all 12001. 
Also, using broad-band excess, iShim et aLl ([20111 ) iden- 
tified strong Ha emitters in the redshift range of 3.8 < 
z < 5.0 that have mostly EW(Hq)> 200 A. A total of 
176 objects satisfy the high-EW criterion, spanning a 
redshift range of 0.35 < z < 2.3, which results in a sur- 
face density of 1 object per square arcmin. We are also 
interested in the redshift evolution of the space density 
of the high-EW galaxies. In order to compare our results 
with the local universe we analyzed the SDSS sample of 
emission- line galaxies at z < 0.1. We used the spectral 
measurements in catalog of the MPA-JHU DR7 releascfEl 
to derive absolute R magnitudes f rom th e SDSS apparent 
r magnitudes using Uester et al.l ([2005T ) conversion solu- 
tions. A passive evolution correction was applied to the 
WISP sample in order to compare the space d ensities at 
a fixed stellar mass (e.g. iSargent et al.l 12007). We ap- 
plied an absolute magnitude cut at —21 < Mr < —18.5 
in the redshift interval 0.05 < z < 0.1. The same abso- 
lute magnitude cut was applied to the high-EW sample. 
In the end, the space density of galaxies with EW > 200 
A is about 5.5 x 10~ 6 Mpc~ 3 in the SDSS, and about 
7 x 10~ 5 Mpc~ 3 in WISP. This represents an evolution 
of more than a factor of 10 in the number density of 
strong emission- line galaxies between z ~ 1.3 and z ~ 0. 

5. CONTRIBUTION OF EMISSION LINES TO BROADBAND 
PHOTOMETRY 

To quantify the contribution of nebular lines to the 
total broadband flux density, we have calculated syn- 
thetic flux densities from our spectra using the transmis- 
sion curves of Jno and Huo filters and the continuum- 
subtracted emission lines only, which were then com- 
pared to the total flux density. The result is shown in 
Figure [4] where we plot the contribution of the line flux 
density to the Jno and Hhq total flux densities as a 
function of the magnitudes in these filters. The nebular 
contribution can represent most of the broadband flux 
density of the galaxy, i.e. more than 50 %, which trans- 
lates into more than 0.75 mag. The distribution has a 
median value of about 25 %, which represents 0.3 mag 

10 http : //www . mpa-garching . mpg . de/SDSS/DR7/ 




1000 
EW„, 



2000 



Fig. 2. — Rest-frame equivalent width distribution for objects 
with EW > 200 A in the WISP Survey. The total number in each 
bin is divided into the [Oin] A5007 line (presented in blue) and the 
Ha line (presented in red). 
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Fig. 3.— The [Om] A5007/H/3 versus [Sn] (A6717+6732)/Ha di- 
agnostic diagram. Six of our high-EW galaxies are represented with 
blue circles with lcr error bars circles. The solid red c urve is the sep- 
aratio n between star- forming galaxies and AGNi of Ke wlev et al.l 
(2001). The red dashed lines show typical model uncertainties of 
± 0.1 dex 

contribution from the emission lines. The presence of 
such strong emission lines in the rest-frame optical spec- 
tra can have several implications on the study of high- 
rcdshift galaxy populations. 

5.1. Contamination of high- redshift galaxy samples 

In order to select very high-redshift galaxies, the Ly- 
man break technique exploits the continuum drop blue- 
ward of Lya caused by the intervening IGM absorp- 
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Fig. 4. — Contribution of the nebular line flux density to the 
total broadband flux density in the J and H bands as a function 
of the total magnitude for objects with EW re jt > 200 A. This 
quantity is calculated by comparing the synthetic flux density of 
the emission lines only to the total flux density in the filter. The 
contribution to the F110W filter flux density is presented in blue 
and the F140W in red. We plot error bars for objects that have a 
nebular contribution higher than 50 %. For object with less than 
50 % contribution, we plot instead an averaged error bars on the 
top-right of the figure. 



tion. By using a set of broad-band filters, one can sam- 
ple the continuum of the galaxy to search for a detec- 
tion in the red band and an absence of signal in the 
bluer bands. This method has been extensively used 
to select z ~ 7—10 candidates us ing, in particular 



the WFC3 IR camera onboard HST (lOesch etal 
Bunker et al. 2010: Yan ct al. 2010: Bouwens ct al 



2010; 



2011 



Lorenzoni et al.l 1201 ID . Other groups have identified 



high-z sources in the WFC3 imaging data set by using 
SED fitting approach to o btain photometric redsh ift dis- 
tributions (jMcLure et al.ll2010t iFinkelstein et aTll2010fl . 
Also, in order to probe the bright end of the luminosity 
function (LF) at z ~ 7, a wide area coverage is neces- 
sary. Two HST observing programs started recently to 
search for bright z > 7 dropout candidates using pure 
parallel WFC3 NIR imaging of ran dom fields: BoRG 
(Brightest of Reionizing Galaxies, iTrenti et al.l 120111) 
and HIPPIES (Hub ble Infrared Pur e Parallel Imaging 
Extragalactic Survey lYan et al1l2011|) . 

However, low-redshift sources can have similar colors to 
those of high-z galaxies and therefore satisfy the dropout 
selection criteria. The possible contamination from cool 
dwarfs, transient objects, old galaxies or spurious de- 
tections have been addressed in previous studies (e.g. 
IBouwens et al.ll2011[ ). We consider in this work another 
potential source of contamination: the strong emission- 
line galaxies at lower redshift. Since the WISP emission- 
line galaxy sample is spectroscopically selected over a 
wavelength range of 0.8 — 1.7/im, we were naturally in- 
terested in the case of z ~> 7 galaxies, so called Y and 
J dropouts. We note that the wide shallow observations 



of ITrenti et all (|201lD and I Yan et all ([20111 ) use the Y 09 $ 
filter, whereas the deep programs use Y105. This can 
slightly change the results presented hereafter because 
of the different width of the filters and the wavelength 
position of the emission lines. 

In order to estimate the contamination of such sources 
to the dropout selection we use the same filter set used 
by the wide surveys, i.e. Y098, J125 and -f/igo, and com- 
pare their colors. The dominating contributions that af- 
fect the broadband colors are [Oiii]AA5007,4959 and Ha 
emission lines (cf. left panel of Figure [SJ . In the right 
panel of Figure [51 we examine the location of the emis- 
sion lines and their relative contribution to the total flux 
density of the filters as a function of redshift. Depend- 
ing upon the wavelength position of the emission lines, a 
galaxy can appear to have a wide range of IR colors and 
scatter into the high-z dropout selection window. 

Figure [5] shows the logs — ^125 vs. J125 — -ffi60 color- 
color diagram commonly used to isolate high-z candi- 
dates populated with the high-EW galaxies. We find 
a large dispersion in both J125 — -H160 an d Y098 — ^125 
colors, and we observe galaxies that fall in the high-z 
selection window (grey region) that denotes the color se- 
lection J 12 f> - #160 < 0.5 and Y 09S - ^125 > 0.5 (e.g. 
IBouwens etalll2010l ). 

The contamination of z ~ 8 galaxy candidates by low- 
z interlopers with stron g emission lines h as be en inves- 
tigated theoretically by iTaniguchi et all (|2010f ). These 
authors used SED models to predict the strength of the 
emission lines from the ionizing continuum, and then ex- 
amined whether their synthetic spectra satisfy the selec- 
tion criterion of z ~ 8 galaxies. The [Oiii] A5007 equiva- 
lent width in their models varies from - 10 3 A to ~ 100 
A for age of 1 Myr to 100 Myr. They conclude that the 
contribution of the emission lines is too small to meet 
the color requirement. We measure EW > 1000 A in 
z ~ 1 — 2 galaxies, however, and these do satisfy the Y- 
dropout selection as we can see in Figure [6] However, one 
of the main conditions used to rule-out low— z interlop- 
ers in such surveys is the non-detection in the optical do- 
main . To this end, HUDF observations (IBouwens et akl 
[20101: iBunker et all [20101 : lYan et al.l I2010D make use of 
ACS data in several optical bands which are at least 1 
magnitude deeper than IR images. This makes the con- 
tamination from low-z strong line emitters very unlikely, 
unless the extinction is very high. We were able to mea- 
sure the reddening from Ha/H/3 ratio for one of the in- 
terlopers of Figure [5] (purple circles) and found a small 
value of E(B-V) - 0.05 (A„ - 0.15). Also, in the rest of 
the high-EW sample, objects for which Ha and H/3 are 
available show little extinction. It is also possible that 
this is the result of a selection effect in the sense that 
we are missing the dusty galaxies where H/3 line remains 
undetectable because it would be fainter than our flux 
limit. 

On the other hand, the wide pure parallel surveys 
BoRG and HIPPIES do not use such deep optical data 
to discard these interlopers. Indeed, the UVIS obser- 
vations in the F606W or F600LP bands reach a depth 
comparable to the IR. The color selection criteria used 
in these studies are Y098 ~ ^125 > 1. 75 and J125 — - ffieo < 
0.02 + 0.15 x (Y 098 - Ji25 -1.75) for ITrenti et all (|20Tl . 
and Y 98 - ^125 > 1-35 and J125 - -ffi60 < 0.3 for 
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Fig. 5. — Contamination of the selection of high-redshift dropout galaxies by lower-redshift sources with strong emission lines. Left: 
Example of a high-EW emission-line object, WISP42.115 at redshift z = 1.36. The throughput curves of the 098A/, F125W, and Fl&OW 
filters used for the synthetic flux calculations are over plotted. Right: Yo98 — J125 vs J125 — #160 color track of strong emission-line galaxy 
as a function of redshift. An example of a high-EW galaxy spectrum is used to calculate synthetic magnitudes in the IR bands, with 
EW([On] A3727, [Om] A5007, Ho) ~ (200, 800, 550). Beyond the wavelength limits of the WISP spectrum, a flat continuum in /„ is 
assumed. 
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Fig. 6. — Yo98 — J125 vs. J125 — -H16O color-color diagram for 
the high-EW sample of galaxies. The magnitude in each band is 
derived from synthetic flux calculation using the grism spectra and 
the WFC3 filter curves. The shaded region denotes the selection 
window of z ~ 8 Y"o98-dropout surveys, which satisfy Yn98 — J125 > 
0.5 and J125 — i?i60 < 0-5- The galaxies marked in blue circles 
satisfy this general color-color selection criteria, whereas the rest 
of sample is shown with black points. In addition, the galaxy 
represented with a bigger blue circle satisfies Yo98 — J125 > 1.3 
and J125 — #160 < 0.3 used by the wide sur veys. We compare this 
galaxy with the Yo gs-dropout candidates of lYan et al.l 1120111) and 
jTrenti et al.ll2l)TlT) , shown in orange and green stars, respectively, 
and whose color criteria are detailed in the text. The error bars 
indicate la uncertainties, while the arrows represent lower limits. 



lYan et all (|2011D , respectively. We plot in Figure [5] the 
colors of the high-EW galaxies (in blue) together with 
the samples of z ~ 8 candidates selected in the two par- 
allel surveys (orange and green). When applying the 
respective color criteria of these two surveys, and given 
the uncertainties, we end up with 1 interloper in both 
the BoRG and HIPPIES samples. In this section, we 
used 24 fields that have both G102 and Gm observa- 
tions, which yields an effective survey area of about 80 
acrmin 2 . When we scale this area to the above surveys, 
we find that the contamination from such sources repre- 
sents about 1 object in every 17 fields observed by the 
wide shallow surveys. This is assuming the same number 
density of such sources at the depth of the wide surveys, 
because their photometric observations go deeper than 
our spectroscopic survey. 

To show the necessity of deep optical imaging, we have 
estimated the expected flux in the optical bands of these 
interlopers using the best fit SEDs of the high-EW galax- 
ies presented in Sect. 15.21 Using the filter through- 
puts, we calculated the magnitudes in the F606W and 
F600LP optical bands. We find that almost all the 
galaxies have red Vgoe — ^098 colors ranging from 0.1 to 
0.9 mag. The same red colors are obtained when we use 
the F600LP filter. These results are also confirmed by 
the g' — J and i' — J colors observed for these objects 
(see Table \2§ ■ While some of the interlopers that sat- 
isfy the IR color criteria can be identified in the UVIS 
filters, a mean value of VI306 — 5"098 ~ 0.6 in our results 
suggests that ~ 1 mag deeper imaging is needed in the 
optical in order to rule out the contamination from all the 
low-z high-EW galaxies. In this sense, the ongoing wide 
survey of the Cosmic Assembly N ear-IR Deep Extra - 
galactic Legacy Survey (CANDELS, IGroein et alJl20ll 
will have appropriately deep optical imaging in multiple 
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HST/ ACS (Advanced Camera for Surveys) bands. 

The high-z narrow band surveys of Lya emitters can 
also be, to some extent, subject to the contamination of 
these strong low-z emitters. The selection method relies 
on the flux excess detected in the narrowband relative 
to the broadband flux, and selects candidates that have 
rest-frame equivalent width higher than ~ 40 A. Ob- 
viously, the galaxies presented in this work have much 
larger EW values, and the narrowband excess can be at- 
tributed in some cases to these lower redshif t interlopers. 
Howev er, a t redshift z ~ 6.5 for in stance, lOuchi et al.1 
(pOTOh and iKashikawa eTaD (pOllT ) use also the color 
blueward of the Lya emission to detect the Lyman break 
of the NB921-selected galaxies. This additional criterion 
significantly reduces the low-z cont amination fr a ction . 
Using R - i > 1 and i - z > 1.2, lOuchi et all (|2008l 
I2010T ) found in their spectroscopic follow-up that such 
contamination is very small for their sample of z ~ 5.7 
and z ~ 6.6 candidates, respectively. Using similar con- 



straints, |Kaihi^^aiFal] (|2011D report 81% (70%) con- 
firmation rate for their z = 6.5 (5.7) LAE candidates. 
Combining the color check blueward of Lya with the 
narrowband selection should be sufficient to identify the 
low-z high-EW sources. 

5.2. Effects on Age and Mass Estimate of High-z 
Galaxies 

While most stellar population synthesis models do not 
include the nebular emission lines, some effort has been 
devoted recently to take into account the impact of such 
contribution on t he SED analysis of high-z galaxies. 
lEvles et al.l (|2007[ ) used the UV star formation rate to 
assess the contribution of optical emission lines to their 
Spitzer IRAC photometry and their e ffect on the SEP 
fitting of z ~ 6 galaxies. More recently. iFinkelstein et all 
(|201lD corrected for the optical emissi on lines in their 
fittin g of two LAEs at z ~ 2.3 (see also [McLindc n et al.l 
120111) . iSchaerer fc de Barrosl (|2009f ) treated the effect of 
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nebular emission (lines and continua) in a self-consistent 
way by predicting the absolute line intensities from the 
ionizing continuum of the template SED. They showed 
that neglecting the nebular component could lead to 
overestimate the age of z ~ 6 galaxies by a factor of 
3. This theoretic al prescription has been followed in 
several pape rs fee lOno et aJ J |2010U Watson et al.ll20Tol; 
Rafter et all I201tt iFinlator et all 120111: iMcLure et al.l 
20111: llnoud 



2010; 

20111 ) dealing with SED modeling. 



Using the sample of very-high EW objects, we empir- 
ically study the impact of nebular emission lines on the 
SED fitting of galaxies in general, and on the age and 
mass estimates in particular. We selected 9 galaxies from 
fields WISP5 and WISP7 that have optical follow up in g 
and i bands obtained with the 5m Hale telescope at Palo- 
mar observatory, in addition to the J and H photometry 
from the WISPS IR observations. 

To model the stellar populati on properties of these 
galaxies, we use th e FAST code (|Kriek et alj |2009) to fit 
iBruzual fc Chariot] (|2003f ) synthesis models to the con- 
tinuum magnitudes. We assume exponentially declining 
star formation history, a metallicity of Z = 0.02, and 



a iChabrierl (|2003l ) initial mass function (IMF). The fit- 
ting results are presented in Figure [TSJ [3 and El The 
black curve is the model spectrum fit to the raw ob- 
served magnitudes. The blue curve is the model fit to 
the magnitudes corrected for the contamination of neb- 
ular emission lines following the procedure detailed in 
Section [5.11 In addition, we used the LRIS spectrum to 
correct the i' band magnitude of WISP5_230. Depend- 
ing on the wavelength location and the strength of the 
contributing emission lines, the difference between the 
models varies dramatically. 

In Table [21 we compare the stellar populations proper- 
ties derived from the two model fits. It is shown that not 
accounting for the contribution of emission lines can sig- 
nificantly affect the age and the stellar mass estimates in 
these galaxies. Specifically, the age and the mass of indi- 
vidual objects can be overestimated by an order of mag- 
nitude. We show for instance that, based on the observed 
magnitudes, the galaxy WISP5_170 has an age of t ~ 160 
Myr and a mass of M* ~ 8 x 10 s M Q . Once corrected 
for the nebular emission, the age and mass are revised to 
t ~ 10 Myr and M* ~ 10 8 Mq, respectively. Considering 
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TABLE 2 

Broadband Photometry of the Sample 



Object z g' i' J Jcorr H H corr 

(mag) (mag) (mag) (mag) (mag) (mag) 

WISP5.53 1.94 22.98 (0.02) 22.52 (0.04) 22.43 (0.14) 22.46 22.14 (0.01) 22.30 

WISP5.73 0.88 24.30 (0.06) 23.28 (0.05) 23.04 (0.19) 23.58 23.03 (0.01) 23.37 

WISP5.108 1.36 24.76 (0.08) 24.11 (0.10) 23.62 (0.25) 23.85 23.57 (0.02) 24.03 

WISP5.170 1.62 25.36 (0.14) 25.19 (0.28) 24.54 (0.38) 24.90 24.61 (0.04) 25.33 

WISP5.217 1.93 24.93 (0.09) 24.55 (0.16) 24.97 (0.45) 25.04 24.12 (0.02) 24.74 

WISP5.230 0.70 25.93 (0.23) 25.84 a (0.16) 25.12 (0.48) 26.17 25.52 (0.08) 25.53 

WISP7.140 1.20 25.00 (0.10) 24.70 (0.18) 23.99 (0.29) 24.52 24.05 (0.04) 24.24 

WISP7.145 2.04 24.65 (0.07) 24.24 (0.12) 24.04 (0.30) 24.24 23.67 (0.02) 23.85 

WISP7.152 1.22 25.03 (0.10) 24.43 (0.14) 24.09 (0.30) 24.46 24.07 (0.03) 24.32 



Note. - The optical and near-infrared magnitudes of our subsample of strong 
emission-line galaxies. The J and H magnitudes are from our WFC3 IR observations, 
while the g' and i 1 data were obtained during our optical follow-up at Palomar obser- 
vatory. The J corr and H corr columns correspond to the magnitude values corrected for 
emission-line contribution. The 1-cr uncertainties are given in parentheses. The ground- 
based photometry aperture was corrected to match that of the WFC3 images (cf. Section 
|3J. Magnitudes are in AB system. 

a The i' band flux of this object was corrected (by 1.3 mag) for the contamination of 
strong emission lines using the LRIS spectrum. 



TABLE 3 

Impact of emission lines on the physical properties 



Object 


■<4 v 


log(SFR) 


log(Age) 


log(Age (corr)) 


log(Mass) 


Mass (corr) 






(Me yr- 1 ) 


(yr) 


(yr) 


(M©) 


(Mq) 


WISP5.53 


1.20 


2.58 


7.00 


7.00 


9.64 


9.52 


WISP5.73 


0.20 


-0.90 


8.20 


8.40 


9.08 


8.91 


WISP5.108 


1.60 


1.83 


7.00 


7.00 


8.89 


8.56 


WISP5.170 


0.70 


0.78 


8.20 


7.00 


8.90 


8.00 


WISP5.217 


0.80 


1.23 


8.00 


7.00 


9.21 


8.35 


WISP5.230 


0.00 


-1.00 


8.30 


8.20 


7.81 


7.67 


WISP7.140 


0.20 


-0.54 


8.10 


7.20 


8.89 


8.43 


WISP7.145 


1.00 


1.63 


7.90 


7.30 


9.45 


9.09 


WISP7.152 


1.50 


1.57 


7.00 


7.00 


8.55 


8.34 



Note. — Physical properties derived from the SED fitting. The best fit for the 
age and and the stellar mass are presented with (corr) and without including the 
nebular emission lines. 



all the objects with varying emission-line intensities, we 
derive a mean age for the sample of t 

mean r ^~' 

100 Myr and 

a mass of M^ etm ~ 1.4 x 10 9 M Q . Correcting for emis- 
sion lines brings these values down to t mean ~ 55 Myr 
and a mass of M^ earl ~ 7.3 x 10 s M Q . This change repre- 
sents an average c orrec tion factor of 2. For comparison, 
iFinkelstein et al] (|2011[ ) derived for their two objects at 
z ~ 2.5 a correction factor of (2, 6) for the mass and a 
factor of (10, 1) for the age. We note that the extinction 
A v can also be slightly overestimated if no correction is 
applied. 

It is clear that not accounting for nebular emission 
lines when modeling their SEDs can introduce large er- 
rors in the derived age and stellar masses of star-forming 
galaxies at all redshifts. Deriving the physical proper- 
ties of higher redshift galaxies is much more problematic. 
Galaxies at z > 6 will be younger, more metal-deficient, 
and will exhibit higher equivalent widths. Moreover, the 
contribution of emission lines to the broadband filters in- 
creases as (1 + z). However, this effect is mitigated by 
the fact that commonly used broad photometric bands 



increase in width when going redder. This will dilute 
the line in this large wavelength band and decrease the 
EW contribution. This empirical demonstration shows 
the need for careful SED modeling by systematically ac- 
counting for nebular emission lines. 

6. HIGH SPECIFIC STAR FORMATION RATE GALAXIES 

The correlation between the SFR of galaxies and 
their stellar mass M* or their luminosity and metallic 
ity has been studied extensively, especially at low red 
shift. Observ ationally, it has been investigated at differ 



ent redshifts jBrinchmann et al 



iNoeske et al 



iDamen et al 



2007b: Elbaz et al 



20091 lLabbe et al 



2004: lErb et alJl2006a: 



200 



3L= 

l IDaddi et al.112007 



20 ID). This SFR— M* 



relation ship is sometimes c alled the galaxy main se- 
quence (jNoeske et al.ll2007b[ ) as a reference to the stellar 
main sequence. While a SFR— M* relation is exhibited 
at all redshifts, it evolves in the sense that galaxies of a 
given stellar mass have higher SFRs at higher redshifts. 
The evolution of the SFR as a function of M* is of par- 
ticular interest as a test of galaxy evolution models. The 
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Fig. 9. — The star formation rate as a function of the stellar 
mass. The SFR is derived from the optical emission lines and 
the stellar masses from the SED fitting of the optical and NIR 
photometry (cf. text for details). Objects located at z > 1.3 are 
shown with red circles and objects at z < 1.3 with blue circles. 
The lcr uncertainties on the mass and SFR are also overplotted. 
The lines denote previous results of several studies of t he SFR-mass 
correl ation: the orange line is the fit to z = 2 sample of Daddi ct al. 
(2007) with the dashed line representing the typical dispersion of 
0.16 d ex. The black line is for galaxies at z = 1 from El baz et al.l 
(2007) with a dispersion of 0.3 dex (dash ed black line) and th e 
green line for galaxies at 0.2 < z < 0.7 from Nocskc ct al. (2007b). 
The thick part of each line represents the correlation in the mass 
range explored by each study and the thin line is an extrapolation 
to lower masses. 

SFR-M* correlation is indeed predicted by cosm ologi- 
cal hydrodynamic si mulations (|Finlator et al .1120 111 ), and 
iBouche et ail (pOlOl ) argue that s mooth gas accretio n is 
responsible for this trend (see also lDave et al.lfeOllD 

However, at high-redshift, the SFR— M* correlation 
has been limited to a relatively high mass range. The 
techniques used for the selection of galaxies include 
mostly the optical o r NIR colors and UV continuum (e.g. 



Steidel et al.l 12004 iDaddi et al.l I2004t IPapoyich et al l 
20061: ISteidel et alj|1996t IShaplev eHd1l2005t lErb et al.l 



2006ai ISchiminovich et alj 120071) . Unfortunately, this 



misses a population of star-forming galaxies with faint 
continua. Photometric redshifts have also been used 
to identify high-z galaxies, but remain unreliable at 
faint magnitudes because of photometric uncertainties 
and emission line contribution. While the most massive 
galaxies are offset to lower specific star formation rates 
(sSFR, the SFR per unit stellar mass) in the SFR- M* 
plane, our sample selection targets precisely galaxies with 
higher sSFRs that would lie in the upper left corner of 
this plane. 

In Figure |H1 we present the star formation rate mea- 
sured from emission lines as a function of the stellar mass 
for 9 of the high-EW sources for which we obtained opti- 
cal follow- up imaging. The SFR is derived from the Ha 
flux using IKennicutti (fl9 98) calibration. No dust correc- 
tion was applied, as we were able to measure Ha and 



H/3 lines in 5 of these galaxies that were consistent with 
no, or little, extinction. We corrected the Ha flux for 
[Nil] contamination using the procedure described in Sec- 
tion For high-redshift galaxies selected by their strong 
[Oill] A5007 line, and where the Ha line is outside the 
G141 window, we also have a detectable [On] A3727 line. 
Likewise, we apply the IKennicutti ()1998[ ) prescription to 
convert the [On] flux into SFR. The stellar masses were 
derived from SED fitting of t he optical and NIR obser- 
vations of the galaxies using IBruzual fc Charlotl J2003) 
spectral synthesis models with exponentially declining 
star formation e~ t//T with t and r as free parameters, a 
IChabrierl ([2001 initial mass function, and metallicity of 
Z = 0.02. 

It can be seen from Figure [§] that these galaxies have 
a higher sSFR than the sequence of normal star-forming 
galaxies at the same redshift derived in previous stud- 
ies. They probe a mass range that has never been ex- 
plored before. We divided the high-EW sample into two 
bins of redshifts below and above z = 1.3 draw n in blue 
and red, respectively. While IDaddi et al.1 (|2007l ) derive a 
tight SFR-M* correlation for z ~ 2 galaxies in GOODS 
(Great Observatories Origins Deep Survey), with a dis- 
persion of 0.16 dex in SFR, the present galaxies lie well 
above the 0.5 dex line. They represent extreme out- 
liers to the relationship at z ~ 2, like the sub- millimeter 
galaxies (SMGs) but orders of magnitude smaller and 
less massive, because we probe stellar masses in the 
range 10 7 5 < M* < 10 9 5 M© yr" 1 , much lower than 
IDaddi et al.l (120071) Com pared to the z ~ 1 GOODS 
sample (|Elbaz et al J 120071 ) . for which the SFR-M* rela- 
tionship has 0.3 dex dispersion, the high-EW sample of 
galaxies are offset by more tha n 1 dex from t he med ian 
value of SFR at a given mass. iNoeske et ah! (|2007aD in- 
terpret z ~ 1 galaxies on the SFR-Af * main sequence as 
the early phase of a star formation history (SFH) that 
smoothly declines for ~ 1 Gyr to z ~ 0. However our 
observed sSFRs denote vigorous star formation episodes 
with very rapid stellar mass build-up, doubling total stel- 
lar mass in < 100 Myr. The observed offset of our galax- 
ies from the SFR-M* main sequence is instead consis- 
tent with strong burst episodes rather than a smoothly 
declining SFH. 

7. EXTREMELY LOW METALLICITY GALAXIES 

It has long been known that lower luminosity galax- 
ies tend also to have lower metallicities according to 
the luminosity -metallicity relationsh i p, although with a 



large scatter (Trcmon ti et al 
ISavaglio et aT] [2005; Er b et al 



l200l ISaTzer et al.l 120051 



l2006a|) . Extremely metal- 



poor galaxies (XMPGs) with oxygen abundance of 
log(Q /H) + 12 < 7.65 a s defined bv IKunth fc Ostlinl 
(2000) are very rare (e.g. iMorales-Luis et al.l 120111 ). In 
the local universe, large EWs galaxies ha ve been sur- 
veyed u sing broad-band c olor selection (jBrown et al.l 
120081) . ICardamone et al.l (|2009l ) found compact star- 
forming galaxies in the SDSS, called "green peas", that 
have very blue r — i colors because of their ver y strong 
[Oiil]A5007 emission line. iKakazu et al.l (|2007l ) showed 
that emission line selection is yet more efficient than 
broad-band selecti on and subsequent followup observa - 
tions (e.g. DEEP2 lDavis et al.ll200l iHovos et al.ll2005l) . 

Many of the faint WISP galaxies have strong [Oin] 
A5007 emission, and often weak [On] A3727 emission, 
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Fig. 10. — LRIS rest-frame spectrum of the XMPG WISP5-230 at z ~ 0.7. We observe a very faint continuum with very strong emission 
lines, such as H/3 or [Om]A5007. The inset shows a close-up of the H7 and [Om] A4363 lines. The detection of the [Oni] A4363 auroral line 
suggests a very low metallicity. The spectrum has been smoothed by 3-pix filter. 



with a m edian ratio of ([O i ii] A5007/[Onl A3727) = 2.5 
(see also iHicks et all 120021: iMaier et all 120061: ILv et al.l 
I2007T ). suggestive of low metal abundances (Log(0/H) < 
8.5)). We present here an example of a high-EW galaxy 
that shows a very low metallicity. 

Amongst 5 high-EW galaxies for which we obtained 
LRIS spectra, only one galaxy was at a suitable redshift 
to observe the optical emission lines necessary for the 
metallicity measurement. The galaxy WISP5_230 proves 
to be an XMPG. It lies at a redshift of z = 0.7 and its 
LRIS spectrum (FigurefT0|) confirms very strong emission 
lines seen in the grism spectra, with no detected contin- 
uum. In particular, the [Oiii] A5007 line has an equiva- 
lent width lower limit of EW resl > 860 A, and the H/3 
line has EW rest > 180 A. We also detect the [Om]A4363 
auroral line, the presence of which always indicates a 
low metallicity, with /([Oin]A4363) = 4.0 ± 0.3 x 10~ 18 
erg s _1 cm~ 2 , which is almost half of the flux of the H7 
line. We used the direct method to derive the metallicity 
using the ratio of [Oiii] A5007,4959 and [Om]A4363 that 
allows us to measure the ele ctron temperature *. U sing 
the NEBULAR package in IRAF (|Shaw k, Dufourlll995[ ). and 
assuming an electron density n e of 100 cm" 3 , we com- 
pute an electron temperature * ~ 2 x 10 4 K. We finally 
derive the oxygen abundance following the equations of 
llzotov et all (20061 ): 

12 + logO+/H+ =log + 5.961 - 0.40 log* 

-0.034* + log(l + 1.35a;), (1) 

12 + log 2+ /H+ = log *4959+A5007 + 6.200 + i^fi 



-0.55 log*- 0.014* (2) 

where x = 10"V*-°- 5 . Since the [On]A3727 line fell on 
a sky line, we first use a fit to local Hii galaxies used in 
iHovos et all (pool , which gives the ([Oiii] A5007/[On] 
A3727) ratio as a function of H/3 equivalent width, where 
the typical uncertainty on the [On] flux is about 50%. 
Then, using the typical ratio of 2.5 observed for the 
WISP galaxies we obtain a [On] flux higher by 50% than 
the previous estimate. We eventually adopt the typi- 
cal [Oiii] /[On] line ratio observed in WISP as the final 
value. In fact, the uncertainties do not strongly depend 
on the precise value of the [On] flux. The flux errors 
were propagated during the calculation. Finally, we ob- 
tain an oxygen abundance of 12+log(0/H)=7.47± 0.11. 
For comparison, the solar metallicity of Z = 0.02 trans- 
lates to 12+log(0/H) - 8.9. 

This object is comparable to the mos t metal-deficient 
galaxies found in the USEL sample (|Hu et al.l 120091 ) 
and close to the l owest metallicity star- forming galax- 
ies known, IZw 18 dThuan & Izotovl [2005l) and SBS0335- 
52 (jlzotov et al.ll2006D . We plot the metallicity and the 
absolute B-band magnitude of WISP5_230 in Figure ITT1 
and compare it with results at different redshifts. The 
luminosity-metallicity relation is an important indicator 
of the chemical enrichment of the inter-stellar medium 
and the mass-loss of galaxies. In the same figure, we 
show the relation established for a large sample of low- 
redshift star-forming galaxies in the SDSS (orange lines, 
iTremonti et"aH 120041). at redshif ts 0.4 < z < 1 (blue 
dashed line. ISavaglio et all 120051) . an d for high - redshif t 
LBGs at z - 2.3 (green triangles, lErb et all l2006al) . 

With an absolute magnitude of M B 16.4, WISP5.230 

is comparable to local XMPGs (magenta stars), includ- 
ing the blue compact galaxy SDSS 0809+1729 found by 
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Fig. 11. — The high-EW emission-line selection is uncovering 
extremely low metallicity objects. We show here the position of 
WISP5-230 in the metallicity-luminosity plane, together with re- 
sults from the literature. The solid orange line represents the me - 
dian of the low-redshift sample in SDSS of Tremonti et al. I|2004l 'l. 
while the dashed and dotted orange lines are the contours that en- 
close 68% and 95% of the sa mple. The blu e dash ed line denotes 
the 2 ~ 0.7 galaxy sample of Savaglio et al. (2005) and the green 
triangles the z ~ 2.3 star- forming galaxies of lErb et al.l (2006a) 
divi ded into six bins of lu minosity. Local XMPGs and GRB hosts 
fromlKcwlcy et al. ( 2007j) are shown in magenta stars an d XMPGs 
of the USEL sample is represented by the black solid line l|Hu et al.l 
120091 ). The position of WISP5.230 is marked with a big red circle. 

IKewlev et all (|200l . 

We observe that, at constant luminosity, it has a metal- 
licity ~ 1 dex lower tha n what would be pr edicted 
from local SDSS galaxies (|Tremonti et all 120041 ) . It is 
better reproduced by the black solid line in Figure fTTI 
that represents the same trend for the USEL sample 
(iHu et al.l |2009T). It is also consistent with relationship 
oflSavaglio et al. (2005) derived for z ~ 0.7 galaxies. 

This object also has a stellar mass of Log(M*) 
= 7.67 ± 0.13 M Q ; the lowest value in the sam- 
ple of 9 galaxies for which we performed SED mod- 
eling (see Section I5.2p . The best- fit SED is shown 
in Figure 1121 The low metallicity and low mass of 
this object are con sistent with the wel l known mass- 
met al licity relation (iTremonti et alll2004l iSavaglio et all 
120051; lErb et al.l I2006at IMannucci et al.l 12009ft . In ad- 
dition, IMannucci et all ([2010ft argue that the mass- 
metallicity relation is simply a projection of a more fun- 
damental relation between mass, gas-phase metallicity 
and SFR. The evolution of the mass-metallicity rela- 
tion up to z = 2.5 would be the result of an evolution 
within this fundamental relation, as higher-SFR galaxies 
are being selected at increasing redshifts. Therefore, for 
low-mass galaxies, the metallicity decreases sharply with 
increasing SFR. However, the high-redshift samples are 
still small and there is not enough overlap between SFRs 
of the low- and high-z samples to confirm this scenario. 
Using the fundamental metallicity equation derived in 



IMannucci et al.l (|2010ft we find that, with SFR~ 1.9 
M Q yr -1 , the galaxy presented here would have an oxy- 
gen abundance of about 12+log(Q/H)~ 7.05. Although 
the m ass-metallicity regime explored in M annucci et all 
(2010) is different, it appears that this galaxy does not 
seem to fit in a scenario of no evolution in the M*-Z- 
SFR fundamental plane. A large sample of galaxies is 
needed to test this scenario, and will be the subject of a 
forthcoming paper (Henry et al. in prep). 

Rest-frame wavelength X [urn] 
0.2 0.4 0.6 0.8 1.0 
24.0 L 'J 
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27.0 L , i , , , . i , , , , i . ^_ 
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Observed Wavelength X [M m ] 

Fig. 12. — Spectral energy distribution fitting of XMPG 
WISP5-230. The black points represent the observed magnitudes 
in the g',i',J and H bands. The blue points are the magnitudes 
corrected for the contribution of nebular emission lines. The J and 
H fluxes are corrected using the WFC3 grism spectra, and the g 1 
flux using the LRIS spectrum. The solid black (blue) curves is the 
best fit SED to the r aw (e mission-line corrected) magnitudes using 
Bruzual & Chariot (20 03) models. The SED fitting procedure is 
detailed in Section 15.21 

8. SUMMARY 

The WISP survey (jAtek et al.ll2010ft offers a unique 
opportunity to search for high-redshift strong emission- 
line galaxies, down to a line flux limit of ~ 5 x 10~ 17 
erg s _1 cm~ 2 , regardless of their continuum brightness, 
thanks to unprecedented IR grisms capabilities of the 
new WFC3 instrument onboard the HST. 

We selected in this work a sample of 176 objects in 54 
fields, with rest-frame EW higher than 200 A, spanning 
a redshift range of 0.35 < z < 2.3. This selection results 
in a surface density of 1 object per square arcmin. 

The presence of such strong emission lines in the spec- 
tra of star-forming galaxies at 0.8 < z < 2.35 has impor- 
tant implications for the study of high-redshift galaxy 
population. We first show that the contribution of neb- 
ular lines to the total broadband flux density can be 
more than 1 magnitude, with a median value at 0.3 
mag. We demonstrate that strong emission lines falling 
in the IR broadband filters can mimic the Yogs — J125 and 
■^125 — Hiqq color criteria used to select z ~ 8 galaxy can- 
didates. While, the presence of such interlopers in deep 
HUDF observations is unlikely, the wide WFC3 pure par- 
allel surveys are prone to such contamination because of 
the lack of deep optical observations. 

High-EW emission lines can significantly affect the 
SED modeling of high-z galaxies and consequently the 
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derived physical properties. We show that, when the 
emission lines create the appearance of a Balmer/4000 A 
break, the age and mass can be overestimated by an aver- 
age factor of 2 and up to a factor of 10. More importantly, 
the situation will be much more problematic at higher 
redshift where we expect higher equivalent widths. Our 
results confirm empirically t he theoretical predictions of 
ISchaerer fc de Barrosl ([20091 ) in a large sample of galaxies 
and call for a careful treatment of nebular lines in SED 
fitting of star-forming galaxies at all redshifts. 

The high-EW sample consists of young, unevolved 
galaxies with high specific star formation rates, which 
appear to be extreme outliers of the main sequence of the 
SFR-M* relationship. The high-EW selection can sam- 
ple metal-deficient galaxies. As an example, we obtained 
rest-frame optical spectroscopy with LRIS at Keck for 5 
of these objects, but only one of them was at the right 
redshift for the metallicity measurement. This galaxy is 
located at z ~ 0.7 and is an extremely metal-poor galaxy. 
It follows the met alli city- luminosity relat ion derived from 
the USEL sample dHu et all 12009ft but is also consis- 
tent with z ~ 0.7 galaxies (jSavaglio et all 12005ft . With 



an oxygen abundance of 12+Log(0/H)= 7.47 ± 0.11, 
it is amongst the lowest metallicities in the samples 
of XMPGs UKunth fc Ostlinl [20001: IKakazu et all 127)071 : 
iKewlev et al.l 12007ft . A larger sample of low-metallicity 
galaxies will be analyzed in a forthcoming paper (Henry 
et al. in prep), but we already demonstrate that we can 
find faint metal-poor galaxies in the high-EW sample to 
be targeted in follow-up spectroscopy. 
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the staff members of the Space Telescope Institute for 
their help with the data reduction. We also wish to rec- 
ognize and acknowledge the very significant cultural role 
and reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are 
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